Gene Regulatory Network Inference and Gene Module Regulating Virulence in Fusarium oxysporum.

In this work, we inferred the gene regulatory network (GRN) of the fungus Fusarium oxysporum by using the regulatory networks of Aspergillus nidulans FGSC A4, Neurospora crassa OR74A, Saccharomyces cerevisiae S288c, and Fusarium graminearum PH-1 as templates for sequence comparisons. Topological properties to infer the role of transcription factors (TFs) and to identify functional modules were calculated in the GRN. From these analyzes, five TFs were identified as hubs, including FOXG_04688 and FOXG_05432, which regulate 2,404 and 1,864 target genes, respectively. In addition, 16 communities were identified in the GRN, where the largest contains 1,923 genes and the smallest contains 227 genes. Finally, the genes associated with virulence were extracted from the GRN and exhaustively analyzed, and we identified a giant module with ten TFs and 273 target genes, where the most highly connected node corresponds to the transcription factor FOXG_05265, homologous to the putative bZip transcription factor CPTF1 of Claviceps purpurea, which is involved in ergotism disease that affects cereal crops and grasses. The results described in this work can be used for the study of gene regulation in this organism and open the possibility to explore putative genes associated with virulence against their host.

The Search for Cryptic L-Rhamnosyltransferases on the Sporothrix schenckii Genome.

The fungal cell wall is an attractive structure to look for new antifungal drug targets and for understanding the host-fungus interaction. Sporothrix schenckii is one of the main causative agents of both human and animal sporotrichosis and currently is the species most studied of the Sporothrix genus. The cell wall of this organism has been previously analyzed, and rhamnoconjugates are signature molecules found on the surface of both mycelia and yeast-like cells. Similar to other reactions where sugars are covalently linked to other sugars, lipids, or proteins, the rhamnosylation process in this organism is expected to involve glycosyltransferases with the ability to transfer rhamnose from a sugar donor to the acceptor molecule, i.e., rhamnosyltransferases. However, no obvious rhamnosyltransferase has thus far been identified within the S. schenckii proteome or genome. Here, using a Hidden Markov Model profile strategy, we found within the S. schenckii genome five putative genes encoding for rhamnosyltransferases. Expression analyses indicated that only two of them, named RHT1 and RHT2, were significantly expressed in yeast-like cells and during interaction with the host. These two genes were heterologously expressed in Escherichia coli, and the purified recombinant proteins showed rhamnosyltransferase activity, dependent on the presence of UDP-rhamnose as a sugar donor. To the best of our knowledge, this is the first report about rhamnosyltransferases in S. schenckii.

Universal Primers for Detection of Novel Plant Capsid-Less Viruses: Papaya Umbra-like Viruses as Example.

For diagnosis of positive-sense single-stranded RNA viruses, primers are usually raised against the sequence encoding capsid proteins, since structural proteins are more conserved. This chapter focuses on the design of primers for a group of novel viruses lacking a capsid, known as papaya Umbra-like viruses (unassigned genus) associated with Papaya Sticky Disease, which represent a threat to papaya production. Based on sequence alignments of a region encoding the RNA-dependent RNA Polymerase, universal primers to detect all the known viruses from four countries are proposed. The Forward universal primer can be used in combination with clade- and subclade-specific primers for rapid virus identification. We walk the reader through downloading sequences from nucleotide databases, doing sequence alignments and phylogenetic tree construction to identify conserved and variable regions as valid primer targets; we also show how to design and analyze the primers.

Auxin perception in Agave is dependent on the species' Auxin Response Factors.

Auxins are one of the most important and studied phytohormones in nature. Auxin signaling and perception take place in the cytosol, where the auxin is sensed. Then, in the nucleus, the auxin response factors (ARF) promote the expression of early-response genes. It is well known that not all plants respond to the same amount and type of auxins and that the response can be very different even among plants of the same species, as we present here. Here we investigate the behavior of ARF in response to various auxins in Agave angustifolia Haw., A. fourcroydes Lem. and A. tequilana Weber var. Azul. By screening the available database of A. tequilana genes, we have identified 32 ARF genes with high sequence identity in the conserved domains, grouped into three main clades. A phylogenetic tree was inferred from alignments of the 32 Agave ARF protein sequences and the evolutionary relationship with other species was analyzed. AteqARF 4, 15, 21, and 29 were selected as a representative diverse sample coming from each of the different subclades that comprise the two main clades of the inferred phylogenetic reconstruction. These ARFs showed differential species-specific expression patterns in the presence of indole-3-acetic acid (IAA) and 2,4-dichlorophenoxyacetic acid (2,4-D). Interestingly, A. angustifolia showed different phenotypes in the presence and absence of auxins. In the absence of auxin, A. angustifolia produces roots, while shoots are developed in the presence of IAA. However, in the presence of 2,4-D, the plant meristem converts into callus. According to our results, it is likely that AteqARF15 participates in this outcome.

Data mining of metagenomes to find novel enzymes: a non-computationally intensive method.

Currently, there is a need of non-computationally-intensive bioinformatics tools to cope with the increase of large datasets produced by Next Generation Sequencing technologies. We present a simple and robust bioinformatics pipeline to search for novel enzymes in metagenomic sequences. The strategy is based on pattern searching using as reference conserved motifs coded as regular expressions. As a case study, we applied this scheme to search for novel proteases S8A in a publicly available metagenome. Briefly, (1) the metagenome was assembled and translated into amino acids; (2) patterns were matched using regular expressions; (3) retrieved sequences were annotated; and (4) diversity analyses were conducted. Following this pipeline, we were able to identify nine sequences containing an S8 catalytic triad, starting from a metagenome containing 9,921,136 Illumina reads. Identity of these nine sequences was confirmed by BLASTp against databases at NCBI and MEROPS. Identities ranged from 62 to 89% to their respective nearest ortholog, which belonged to phyla Proteobacteria, Actinobacteria, Planctomycetes, Bacterioidetes, and Cyanobacteria, consistent with the most abundant phyla reported for this metagenome. All these results support the idea that they all are novel S8 sequences and strongly suggest that our methodology is robust and suitable to detect novel enzymes.

A Lineage of Begomoviruses Encode Rep and AC4 Proteins of Enigmatic Ancestry: Hints on the Evolution of Geminiviruses in the New World.

The begomoviruses (BGVs) are plant pathogens that evolved in the Old World during the Cretaceous and arrived to the New World (NW) in the Cenozoic era. A subgroup of NW BGVs, the "Squash leaf curl virus (SLCV) lineage" (S-Lin), includes viruses with unique characteristics. To get clues on the evolutionary origin of this lineage, a search for divergent members was undertaken. Four novel BGVs were characterized, including one that is basal to the group. Comparative analyses led to discover a ~670 bp genome module that is nearly exclusive of this lineage, encompassing the replication origin, the AC4 gene, and 480 bp of the Rep gene. A similar DNA module was found in two curtoviruses, hence suggesting that the S-Lin ancestor acquired its distinctive genomic segment by recombination with a curtovirus. This hypothesis was definitely disproved by an in-depth sequence analysis. The search for homologs of S-Lin Rep uncover the common origin of Rep proteins encoded by diverse Geminiviridae genera and viral "fossils" integrated at plant genomes. In contrast, no homolog of S-Lin Rep was found in public databases. Consequently, it was concluded that the SLCV clade ancestor evolved by a recombination event between a primitive NW BGV and a virus from a hitherto unknown lineage.

Molecular characterization of laccase genes from the basidiomycete Trametes hirsuta Bm-2 and analysis of the 5' untranslated region (5'UTR).

The aim of this study was to identify and characterize laccase genes produced by Trametes hirsuta Bm-2 in a liquid medium, both with and without induction. The amplification of 5'and 3'regions of laccase sequences was obtained by the RACE-PCR method, and these were assembled to obtain a cDNA of total length. Two new laccase genes were isolated from basal medium (lac-B) and lignocellulosic grapefruit substrate (lac-T), both encoding open reading frames of 2566 bp. Both laccase-predicted proteins consisted of 521 amino acids, four copper-binding regions, a signal peptide, and five potential glycosilation sites (Asn-Xaa-Ser/Tre). Moreover, the deduced amino acid sequences share about 76-85% identity with other laccases of WRF. Sequence comparison showed 47 synonymous point mutations between lac-B and lac-T. In addition, 5' untranslated regions (UTR) of laccase genes lac-B and lac-T showed differences in length and number of regulatory elements that may affect transcriptional or translational expression of these genes.

Ploidy Determination in the Pathogenic Fungus Sporothrix spp.

The pathogenic clade of the Sporothrix genus comprises the etiological agents of sporotrichosis, a worldwide emergent disease. Despite the growing understanding of their successful pathogen traits, there is little information on genome sizes and ploidy within the genus. Therefore, in this work, we evaluated the ploidy of four species of the Sporothrix genus, specifically Sporothrix brasiliensis, Sporothrix schenckii, Sporothrix globosa, and Sporothrix pallida. Through cell cycle analysis of the yeast-phase cells, we showed that the DNA content of G0/G1 cells was similar to the genome size determined by whole genome sequencing. Moreover, ploidy of S. schenckii, S. brasiliensis, and S. pallida that was determined by allele composition using next-generation sequencing (NGS) data is consistent with monomorphic positions at each allele. These data show that the analyzed strains of Sporothrix are haploid, or at least aneuploid, thereby laying the foundation for the development of a molecular toolbox for Sporothrix spp.

Characterization of N-Acetylglucosamine Biosynthesis in Pneumocystis species. A New Potential Target for Therapy.

N-acetylglucosamine (GlcNAc) serves as an essential structural sugar on the cell surface of organisms. For example, GlcNAc is a major component of bacterial peptidoglycan, it is an important building block of fungal cell walls, including a major constituent of chitin and mannoproteins, and it is also required for extracellular matrix generation by animal cells. Herein, we provide evidence for a uridine diphospho (UDP)-GlcNAc pathway in Pneumocystis species. Using an in silico search of the Pneumocystis jirovecii and P. murina (Pm) genomic databases, we determined the presence of at least four proteins implicated in the Saccharomyces cerevisiae UDP-GlcNAc biosynthetic pathway. These genes, termed GFA1, GNA1, AGM1, and UDP-GlcNAc pyrophosphorylase (UAP1), were either confirmed to be present in the Pneumocystis genomes by PCR, or, in the case of Pm uap1 (Pmuap1), functionally confirmed by direct enzymatic activity assay. Expression analysis using quantitative PCR of Pneumocystis pneumonia in mice demonstrated abundant expression of the Pm uap1 transcript. A GlcNAc-binding recombinant protein and a novel GlcNAc-binding immune detection method both verified the presence of GlcNAc in P. carinii (Pc) lysates. Studies of Pc cell wall fractions using high-performance gas chromatography/mass spectrometry documented the presence of GlcNAc glycosyl residues. Pc was shown to synthesize GlcNAc in vitro. The competitive UDP-GlcNAc substrate synthetic inhibitor, nikkomycin Z, suppressed incorporation of GlcNAc by Pc preparations. Finally, treatment of rats with Pneumocystis pneumonia using nikkomycin Z significantly reduced organism burdens. Taken together, these data support an important role for GlcNAc generation in the cell surface of Pneumocystis organisms.

The sexual state of Aspergillus parasiticus.

The sexual state of Aspergillus parasiticus, a potent aflatoxin-producing fungus within section Flavi, is described. The production of nonostiolate ascocarps surrounded by a separate peridium within the stroma places the teleomorph in genus Petromyces. Petromyces parasiticus differs from P. alliaceus by its larger ascospores with finely tuberculate ornamentation. The anamorphic Aspergillus states of the two species differ in conidial head color and microscopic characters.

Sexual reproduction and recombination in the aflatoxin-producing fungus Aspergillus parasiticus.

The fungal phylum Ascomycota comprises a large proportion of species with no known sexual stage, despite high genetic variability in field populations. One such asexual species, Aspergillus parasiticus, is a potent producer of carcinogenic and hepatotoxic aflatoxins, polyketide-derived secondary metabolites that contaminate a wide variety of agricultural crops. In this study, individuals of A. parasiticus from a population showing an evolutionary history of recombination were examined for sexual reproduction. Crosses between strains with opposite mating-type genes MAT1-1 and MAT1-2 resulted in the development of ascospore-bearing ascocarps embedded within stromata. Sexually compatible strains belonged to different vegetative compatibility groups. Recombination through the independent assortment of chromosomes 3 and 6 was detected using loci for mating type, aflatoxin gene cluster, and a protein-encoding gene. Our discovery of the sexual stage in A. parasiticus has important implications for current biological control strategies using nontoxigenic strains to reduce aflatoxin contamination in crops.

Characterization and population analysis of the mating-type genes in Aspergillus flavus and Aspergillus parasiticus.

We characterize the mating-type genes in Aspergillus flavus,Aspergillus parasiticus and Petromyces alliaceus. A single MAT1-1 or MAT1-2 gene was detected in the genomes of A. flavus and A. parasiticus, which is consistent with a potential heterothallic organization of MAT genes in these species. In contrast, the only known, functionally homothallic species in Aspergillus section Flavi, P. alliaceus, has tightly linked (<2kb) MAT1-1 and MAT1-2 genes, typical of other self-fertile homothallic euascomycetes. This is the first example of linked MAT genes within a homothallic species of Aspergillus. We tested the null hypothesis of no significant difference in the frequency of MAT1-1 and MAT1-2 in A. flavus and A. parasiticus sampled from a single peanut field in Georgia. For each species, mating-type frequencies were determined for the total population samples and for samples that were clone-corrected based on vegetative compatibility groups (VCGs) and aflatoxin gene cluster haplotypes. There was no significant difference in the frequency of the two mating types for A. flavus and A. parasiticus in either VCG or haplotype clone-corrected samples. The existence of both mating-type genes in equal proportions in A. flavus and A. parasiticus populations, coupled with their expression at the mRNA level and the high amino acid sequence identity of MAT1-1 (77%) and MAT1-2 (83%) with corresponding homologs in P. alliaceus, indicates the potential functionality of these genes and the possible existence of a sexual state in these agriculturally important species.

Recombination, balancing selection and adaptive evolution in the aflatoxin gene cluster of Aspergillus parasiticus.

Aflatoxins are toxic and carcinogenic polyketides produced by several Aspergillus species that are known to contaminate agricultural commodities, posing a serious threat to animal and human health. Aflatoxin (AF) biosynthesis is almost fully characterized and involves the coordinated expression of approximately 25 genes clustered in a 70-kb DNA region. Aspergillus parasiticus is an economically important and common agent of AF contamination. Naturally occurring nonaflatoxigenic strains of A. parasiticus are rarely found and generally produce O-methylsterigmatocystin (OMST), the immediate precursor of AF. To elucidate the evolutionary forces acting to retain AF and OMST pathway extrolites (chemotypes), we sequenced 21 intergenic regions spanning the entire cluster in 24 A. parasiticus isolates chosen to represent the genetic diversity within a single Georgia field population. Linkage disequilibrium analyses revealed five distinct recombination blocks in the A. parasiticus cluster. Phylogenetic network analyses showed a history of recombination between chemotype-specific haplotypes, as well as evidence of contemporary recombination. We performed coalescent simulations of variation in recombination blocks and found an approximately twofold deeper coalescence for cluster genealogies compared to noncluster genealogies, our internal standard of neutral evolution. Significantly deeper cluster genealogies are indicative of balancing selection in the AF cluster of A. parasiticus and are further corroborated by the existence of trans-species polymorphisms and common haplotypes in the cluster for several closely related species. Estimates of Ka/Ks for representative cluster genes provide evidence of selection for OMST and AF chemotypes, and indicate a possible role of chemotypes in ecological adaptation and speciation.

Gene duplication, modularity and adaptation in the evolution of the aflatoxin gene cluster.

BACKGROUND: The biosynthesis of aflatoxin (AF) involves over 20 enzymatic reactions in a complex polyketide pathway that converts acetate and malonate to the intermediates sterigmatocystin (ST) and O-methylsterigmatocystin (OMST), the respective penultimate and ultimate precursors of AF. Although these precursors are chemically and structurally very similar, their accumulation differs at the species level for Aspergilli. Notable examples are A. nidulans that synthesizes only ST, A. flavus that makes predominantly AF, and A. parasiticus that generally produces either AF or OMST. Whether these differences are important in the evolutionary/ecological processes of species adaptation and diversification is unknown. Equally unknown are the specific genomic mechanisms responsible for ordering and clustering of genes in the AF pathway of Aspergillus. RESULTS: To elucidate the mechanisms that have driven formation of these clusters, we performed systematic searches of aflatoxin cluster homologs across five Aspergillus genomes. We found a high level of gene duplication and identified seven modules consisting of highly correlated gene pairs (aflA/aflB, aflR/aflS, aflX/aflY, aflF/aflE, aflT/aflQ, aflC/aflW, and aflG/aflL). With the exception of A. nomius, contrasts of mean Ka/Ks values across all cluster genes showed significant differences in selective pressure between section Flavi and non-section Flavi species. A. nomius mean Ka/Ks values were more similar to partial clusters in A. fumigatus and A. terreus. Overall, mean Ka/Ks values were significantly higher for section Flavi than for non-section Flavi species. CONCLUSION: Our results implicate several genomic mechanisms in the evolution of ST, OMST and AF cluster genes. Gene modules may arise from duplications of a single gene, whereby the function of the pre-duplication gene is retained in the copy (aflF/aflE) or the copies may partition the ancestral function (aflA/aflB). In some gene modules, the duplicated copy may simply augment/supplement a specific pathway function (aflR/aflS and aflX/aflY) or the duplicated copy may evolve a completely new function (aflT/aflQ and aflC/aflW). Gene modules that are contiguous in one species and noncontiguous in others point to possible rearrangements of cluster genes in the evolution of these species. Significantly higher mean Ka/Ks values in section Flavi compared to non-section Flavi species indicate increased positive selection acting in the evolution of genes in OMST and AF gene clusters.